Nonlinear interaction of charged particles with a free electron gas beyond the 

random-phase approximation 
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A nonlinear description of the interaction of charged particles penetrating a solid has become of 
basic importance in the interpretation of a variety of physical phenomena. Here we develop a many- 
body theoretical approach to the quadratic decay rate, energy loss, and wake potential of charged 
particles moving in an interacting free electron gas. Explicit expressions for these quantities are 
obtained either within the random-phase approximation (RPA) or with full inclusion of short-range 
exchange and correlation effects. The Z\ correction to the energy loss of ions is evaluated beyond 
RPA, in the limit of low velocities. 
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When charged particles pass through a solid, energy 
can be lost to the medium through wirious types of elas- 
tic and inelastic collision processes.^ While at relativis- 
tic velocities radiative losses may become important, for 
moving charged particles in the non-relativistic regime 
the energy loss is primarily due to electron-electron (e-e) 
interactions giving rise to the generation of electron-hole 
pairs, collective excitations such as plasmons, and inner- 
shell excitations and ionizations. Energy losses due to 
nuclear recoil are negligible, unless the projectile veloc- 
ity is very small compared to the mean speed of electrons 
in the solid. 

The inelastic decay rate of charged particles in a de- 
generate interacting free electron gas (FEG) has been 
calculated for many years in the first-Born approxima- 
tion or, equivalently, within linear-response theory. This 
is a good approximation when the velocity of the pro- 
jectile is much greater than the average velocity of tar- 
get electrons. However, in the case of projectiles moving 
with smaller velocities, nonlinearities have been shown to 
play a key role in the interpretation of a variety of exper- 
iments. Energy-loss measurements have revealed differ- 
ences, not present within linear-response theory-Jpetween 
the energy loss of protons and antiprotonsBa More- 
over, experimentally observed coherent double-plasmon 
excitationsoQ cannot be described within linear-response 
theory, and nonlinearities may also play an important 
role in .the electronic wake generated by moving ions in 
a FEG.Q 

Pioneering nonlinear calculations of the electronic en- 
ergy loss of low-energy ions, in an electron gas were per- 
formed by Echenique et al u These authors computed the 
scattering cross section for a statically screened potential, 
which was determined splf-consistently using density- 
functional theory (DFT) These static-screening calcula- 
tions have recently beep-extended to velocities approach- 
ing the Fermi velocity!!!] Second-order perturbative cal- 



culations, which do not have the limitation of being re- 
stricted to low projectile velocities, have been reported by 
different authors with use of the random-phase approxi- 
mation (RPA) and by treating the moving charged I?sf U-i 
cle as a prescribed source of energy and momentum.tlrfij 

In this paper, we report a many-body theoretical ap- 
proach to the quadratic decay rate, energy loss, and wake 
potential of charged particles moving in an interacting 
FEG. The decay rate is derived from the knowledge of the 
projectile self-energy. The energy loss and wake poten- 
tial are obtained within quadratic response theory. While 
the first-order contribution to the energy loss may also be 
obtained from the imaginary part of the projectile self- 
energy by simply inserting the energy transfer inside the 
integrand of this quantity, our results indicate that this 
procedure cannot be generalized to the description of the 
second-order energy loss. Unless otherwise is stated, we 
use atomic units throughout, i.e., e 2 = % = m e = 1. 

We consider the interaction of a moving probe parti- 
cle of charge Z\ and mass M with a FEG of density n. 
The probe particle is assumed to be distinguishable from 
the electrons in the Fermi gas, which is described by an 
isotropic homogeneous assembly of interacting electrons 
immersed in a uniform background of positive charge and 
volume V . Many-body theory shows that the probability 
for the probe particle to occupy a given excited state of 
four-momenturp_p decays exponentially in time with the 
decay constantly 

r -1 = -2Im£ P! (1) 

where S p represents the particle self-energy. 

It is well known that the self-energy cannot be com- 
puted by simply evaluating the lowest-order Feynman di- 
agrams, because of severe infrared divergences due to the 
long-range Coulomb interaction. Instead, one needs to 
resum electron-loop corrections and expand in terms of 
the dynamically screened interaction. Up to third order 
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in Z\ , the self-energy of the probe particle can be repre- 
sented diagrammatically as shown in Fig. 1. The sum 
of the first two diagrams represents the so-called GW 
approximation, and the third diagram accounts for Z\ 
corrections. One finds 
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(7) 



The decay rate of Eq. (|6|) has not been reported before. 

In the RPA, density correlation functions are obtained 
by summing over all ring-like diagrams, 



/dq 4 
7^3J U q D p - q [(1+Xq «q) 



RPA 

Aq 



4 + X° q v^xf PA 



-2iZi 



(27T) 



4 Dp—q! Dp— q j r q x ig,— U qi W q _ qi 



and 



, (2) 



Y 



RPA 



K, 



RPA yO tsRPA ts-RPA 
91 91.92 -92 -93 ' 



(8) 



(9) 



where u q is the Fourier transform of the bare Coulomb 
potential, D p is the probe-particle propagator 
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DP P ~ W p - Ej; 
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where x q an d Y qi q2 represent noninteracting density 

correlation functions. Improvements on the RPA are 

typically carried out by introducing an effective e-e 
interaction!!!! 



uj p = p 2 /(2M) is the noninteracting energy, and q is a 
positive infinitesimal. Xq an d Y qi m represent exact time- 
ordered density correlation functions of the interacting 
FEG, 
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where G q is the so-called local-field factor accounting 
for short-range exchange and correlation (xc) effects not 
present in the RPA. Accordingly, the density correlation 
functions Xq an d Y qi q2 are found to be of the RPA form, 
but with all e-e bare-Coulomb interactions Vq replaced 
by v q , i.e., 
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^ the test_charge-electron inverse 
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(pq) nl being the matrix element of the Fourier transform 
of the electron-density operator, taken between exact 
many-electron states of energy E n and E\ , uj n i = E n — Ei , 
and q 3 = -{q x + q 2 ). 

If the probe particle is an ion (M >> 1), the propaga- 
tor D p and the energy p° entering Eq. (|^) can be safely 
approximated by the noninteracting propagator D® and 
energy u> p . Recoil can also be neglected, and the intro- 
duction of Eq. (||) into Eq. (|TJ) then yields, after some 
work of rearrangement, the following expression: 

T-'=A«Zlj |jj, q % 8 -q.v)e( ? ») 
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In the RPA, K q and K q coincide. 

The potential induced in a uniform FEG by the pres- 
ence of the recoiless probe particle may be obtained with 
the use of time-dependent perturbation theory. Keeping 
terms of first and second order in the external perturba- 
tion, one finds 

V ind (r,t) = 2nZ 1 J ^ e < (""-« D *) Wq 
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where v is the particle velocity, <d(x) represents the Heav- 
iside step function, and K q is the so-called test_charge- 
test.charge inverse dielectric function: 



xu qi u q _ qi 5 (g? -qi -v)] , (14) 

where we have introduced the retarded counterparts of 
the time-ordered functions K q and Y q _ 9l entering Eq. 
®. 

The average energy lost per unit length traveled by the 
probe particle is obtained as the retarding force due to 
the potential of Eq. (|l4|) induced in the vicinity of the 
projectile itself. One easily finds 
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Both Eq. (||) for the decay rate and Eq. ( |l5| ) for the 
energy loss can be derived, within many-body perturba- 
tion theory, from the knowledge of the probability for 
the probe particle to transfer four-momentum q to the 
FEG.E3 Nevertheless, this probability cannot be identi- 
fied with the integrand of Eq. @, and therefore the 
energy loss of Eq. (|T^) cannot be obtained by simply 
inserting the energy transfer inside the integrand of Eq. 
(§). 

In the framewor-k_of time-dependent density-functional 
theory (TDDFT),B quadratic response theory yields 
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where an d ^gi,'g 2 represent retarded noninteracting 
density correlation functions, and denotes the Fourier 
transform of 
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V" a:c ([n], x) ibeing the exact time-dependent xc potential 
of TDDFTM Within RPA f* c = 0, and introduction of 
Eqs. @ and (||) into Eqs. (O) and (jl|) then yields 



the results derived in Refs. 
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At this point, we present an application of our formal- 
ism, namely the low- velocity limit of the energy loss with 
inclusion of short-range xc effects. For low projectile ve- 
locities (v -> 0), only the static (ui -> 0) Xq'°, Y m'^ 
and c enter in the evaluation of the energy loss of Eq. 
( |l5| ) , which is then easily found to be proportional to the 
projectile velocity. In particular, in the so-called local- 
density approximation (LDA), which is rigorous in the 
long- wavelength limit (q — * 0), one finds 
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E c (n) being the correlation energy of a uniform electron 
gas of density n. 

Diffusion Monte Carlo calculations of y9 have shown 
that the LDA static xc kernel of Eq. ( plf ) -reproduces 
correctly the static response for all q < 2qpSB We have 
calculated static density correlation functions from Eqs. 



(U) and (0), with use of the LDA static xc kernel of 
Eq. (|2l| ) and the Perdew-Zungem parametrization of 
the quantum Monte—Carlo correlation energy E c (n) of 
Ceperley and Alder .cfl 

In Fig. 2 we show linear (cx Z 2 ) and quadratic (cx Zf) 
contributions to the low-velocity energy loss of Eq. (|15| ) 
divided by the velocity of the projectile, as a func- 
tion of the electron-density parameter r s EB Comparison 
between our new results (solid and short-dashed lines) 
and those obtained within RPA (long-dashed and dotted 
lines) indicates that xc effects become increasingly im- 
portant as the electron density decreases, the impact of 
these effects being more pronounced for the Z\ than for 
the Z\ contribution to the energy loss. As a result, the 
importance of the quadratic contribution increases when 
xc effects are included, and for r s ~ 2.5 it is equal in 
magnitude to the linear contribution [within RPA, linear 
and quadratic contributions are equal for r s ~ 5] . 

The crosses and rhombs in Fig. 2 represent the full 
nonlinear contribution to the energy loss [difference be- 
tween the total energy loss and the linear contribution], 
multiplied by a factor of —1, as obtained frmn DFT cal- 
culations for antiprotons to all orders in Z\ p3 with xc ef- 
fects excluded (crosses) and with LDA xc effects included 
(rhombs). Both full nonlinear calculations are found to 
agree nicely with our quadratic-response calculations in 
the high-density limit (r s — > 0). As the electron density 
decreases, quadratic-response calculations overestimate 
the energy loss of antiprotons, especially when xc effects 
are included, showing that xc and nonlinear (beyond Z\ ) 
effects tend to compensate. 

Finally, we note that linear and quadratic contribu- 
tions to the energy loss can be extracted from the small- 
Z\ behaviour of full nonlinear DFT calculations, which 
have been reported in the low- velocity limit .El When 
(—dE/dx)/Z 2 is plotted as a function of Zi, the re- 
sult has a linear Z\ dependence at small Z\. The in- 
tercept and the slope of this curve at Z\ = give lin- 
ear and quadratic contributions to the full nonlinear en- 
ergy loss, as shown in Ref. [l^ by simply ignoring xc ef- 
fects. These authors repeated their full nonlinear DFT 
calculations with LDA xc included, but were unable to 
compare them with quadratic-response calculations that 
included xc effects at the same level of approximation. 
The results of these quadratic-response calculations are 
now shown in Fig. 2. It can be seen that the ratio 
between xc-included (short-dashed line) and xc-excluded 
(dotted line) quadratic contributions to the energy loss 
is in excellent agreement with the ratios reported in Ref. 
[l2] [see Fig. 7 of this reference] as derived from the in- 
dependence of full nonlinear DFT calculations. 

In conclusion, we have developed a many-body theo- 
retical approach to the quadratic decay rate, energy loss, 
and wake potential of charged particles moving in an elec- 
trons gas, with full inclusion of short-range xc effects. 
We have shown that in the limit of high electron densi- 
ties and low projectile charges our calculated quadratic 
energy loss, which can be extended to the case of larger 
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velocities, reproduces DFT calculations for antiprotons, 
as long as exchange and correlation are treated at the 
same level of accuracy. 
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FIG. 1. The probe-particle self-energy, up to third order 
in Zi. Thick solid lines represent the exact probe-particle 
propagator, iD p . Dashed lines represent the bare Coulomb 
interaction, — iv^. Two- and three-point loops represent 



time-ordered density correlation functions, i\q and 
respectively. 
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FIG. 2. Z\ and Z\ contributions to the low-velocity elec- 
tronic stopping power of Eq. ( |l5| ) divided by the velocity of 
the projectile, as a function of r s , for Z\ = 1. Our full cal- 
culations, as obtained with inclusion of LDA exchange and 
correlation, are represented by solid (Zf contribution) and 
short-dashed {Z\ contribution) lines. Long-dashed {Zi) and 
dotted (Z\ ) lines represent the corresponding results obtained 
within RPA (fq C = 0). The difference [multiplied by a factor 
of —1] between the full nonlinear stopping power for antipro- 
tons, as reported in Ref. ^6| and the linear contribution is 
represented by crosses (xc effects excluded) and rhombs (LDA 
xc included). 
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